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Abstract 

We investigate the twist-3 pion and kaon distribution amplitudes of the pseudoscalar (0^ ^) and 
\^ ' pseudotensor {(j)^ ^) types, based on the effective chiral action from the instanton vacuum. Flavor 

' SU(3) symmetry breaking effects are explicitly taken into account. The Gegenbauer moments and 

the moments of the distribution amplitudes ((C™)) are also computed. Our results are summarized 
as follows: ~ 0.4, af ~ 0.02 and a2i^ ~ 0.14, and ~ 0.02 and ~ ~ 0; 

O ■ (^2^P ~ {CYk^ 0-37 and'(0^ ~ 0, and (C^)^ ~ (C^)^ ~ 0.20 and (^^ ~ 0. We compare our 

results with those from the QCD sum rules. We also discuss the relevant Wilson coefficients which 
were analyzed recently in chiral perturbation theory. 
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I. INTRODUCTION 



The meson light-cone distribution amphtude (DA) plays an essential role in describing 
exclusive hadronic reactions While the leading-twist DAs provide major contri- 

bution to those processes, higher-twist DAs are suppressed by the inverse of the momentum 
transfer Q^, so that higher-twist DAs are less significant in studying high-energy exclusive 
processes. Because of this fact, not many investigations on higher-twist DAs were per- 
formed 0, S H, S EO; III, compared to the leading- twist ones. However, higher- twist 



DAs should be still important, since they give crucial information on the soft part of ex- 
clusive hadron reactions, in particular, in the smaller region. Moreover, it describes the 
transverse motion of quarks inside the meson 0, ^| . 

In the present work, we aim at investigating the two-particle twist-3 pion and kaon 
DAs within the framework of the nonlocal chiral quark model (xQM) from the instanton 
vacuum, taking into account SU(3) symmetry breaking effects. The instanton model of 
the QCD vacuum provides a good framework to study the structure of light mesons, since 
spontaneous chiral symmetry breaking is well realized in the instanton vacuum via the quark 
zero mode ^ 3, ^| ■ This instanton vacuum model was later extended by introducing the 
current-quark masses [l^ IE 13 ■ It was assumed in the model that the large Nc expansion 
is the reasonable one and the results were obtained in the leading order in this expansion. 
In the present approach, we employ the modified effective chiral action with flavor SU(3) 
symmetry breaking effects derived from the instanton vacuum 3, IE 17 L_ 



The model has been applied to describe the leading-twist meson DAsiSlIiHIiilil. 
One of the two-particle twist-3 pion DAs was already studied within the nonlocal xQM in 



the chiral limit [2l| . There are two independent particle-two twist-3 distribution amplitudes 
defined as: 
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X {0\,ljf{Th)t{pfi - P ■ nh,^g{-rn)\M{P)), (2) 



where the subscript Ai denotes either the pion {Ai = vr) or the kaon [Ad = K) with 
on-mass-shell momentum P in the light-cone frame, u and n'^ stand for the longitudinal 
momentum fraction and the light-like vector satisfying = 0, respectively. The spatial 
separation between the quarks inside the meson is represented by r ■ n. In the present work, 
we drop the Wilson line by considering the light-cone gauge, A ■ n = 0. mj and nig are the 
corresponding current quark masses with given flavors / and g. As for the present case, we 
assign them as {f,g,Ai) = {s,u,K^) and (rf, M,7r+). m_vi denotes the corresponding meson 
mass. -Fa4 is the corresponding meson decay constant which is the normalization constant 
for the DAs. The DAs of Eqs. (Q) and (j2I) satisfy the following normalization conditions: 

du(j)%{u)= du(j)'^{u) = 1. (3) 

^0 

The present work is organized as follows: Section II is devoted to the general formalism. 
The numerical results and discussions are given in Section III. We summarize and draw 
conclusions in Section IV. 
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II. GENERAL FORMALISM 



In order to evaluate the nonlocal hadronic matrix elements of Eqs. (0) and (0), we start 
from the low-energy effective QCD partition function derived from the instanton vacuum in 
Euclidean space as follows BMEQ: 



Z 



iil)'^j{x){i(f) + imf) ipf{x 

d^k d^p 
(27r)8 



Mf{k)M,{p)ijf{k)UlUp) 



(4) 



where ip and M. are the fields for the quark and pseudoscalar meson. Mf{k) is the dynam- 
ically generated quark mass discussed later. The background pseudoscalar meson field If^^ 
is given by 



75 



1 + — 75A<"A" - 



1 



a\2 



(5) 



A"^ is the well-known SU(3) Gell-Mann matrices. Now we are in a position to construct 
the quark propa gato r in the instanton ensemble. There have been several approaches for 
this purpose 



gato r m 

15L IigI. [itI . In the present work, we make use of the prescription suggested 
23[ and Musakhanov [J^]. We attempt to discuss it briefly in order. First, 



by Pobylitsa 

Diakonov et al. introduced the quark propagator under the instanton effects by taking into 
account the following assumption 
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ipi{x)il)\{x 
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where Aj and iplx) stand for the single instanton field and the quark zero mode, respectively. 
However, this assumption loses information on higher quark loops although it turns out that 
it works very well from phenomenological point of view. On the contrary, Pobylitsa expand 
the quark propagator {x\{if) + /!/ + im)~^\x) straightforwardly by virtue of the large Nc 
limit [23|. Thus, all non-planar diagrams are removed from the expansion. One then obtains 
an integral equation as follows: 



{x\ 



+ /ij + im 
1 



\x) 



{x 



if + imf 



x) H tr^ 

' 2VN, ' 



d^Zi 



{x\{i^ + J\i + imf)\x) - - 



-1 



(7) 



where j indicates the instanton coordinate for the instanton and antiinstanton. By solving 
Eq. d?!) using an ansatz a{k) for the full propagator, [id + inu + a{id)]~^, one can obtain 
the current-quark mass dpependent dynamical quark mass [la, |2^ . Final expression for it 
can be written as follows: 



Mf{k) = MoF^kp) 



m 



1 + ^-^ 



d^ 



d 



where d 



0.08385 Svrp 

2Nr i?2 



0.198 GeV. 



(8) 
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We notice that, as a consequence, all the effects from the current-quark mass are included 
in the bracket of Eq. (jH}. Here, we employ the standard values for the instanton ensem- 
ble, N/V = 200^ MeV^ and p ~ 1/3 fm ~ 1/600 MeV"^ for the numerical calculations. 
From these values, we obtain Mq = 0.350 GeV. Although the form factor F{k) can be de- 
rived from the Fourier transform of the quark zero mode, we will use the following simple 
parameterization for it: 



Mf{k) = Mo 



(nA2 - A;2 + ie) 



2n 



1 + 



m. 



d? 



d 



(9) 



Note that we are now working in Minkowski space, in which the DA are well defined, rather 
than in Euclidean one. In the present work as done in our previous work [i^, we assume 
that the partition function of Eq. can be analytically continued to Minkowski space. 
Although we do not prrovide a firm theoretical proof for this assumpmtion, it turns out 

A is 
We use 



21 



that it works qualitatively well from the phenomenological point of view [19|, |2 
the cutoff mass which can be taken as the scale par ameter of the present work. 
A ~ 1.0 GeV for the numerical calculation as in j22|. In order to check the effect of the 
parameterization of the form factor, we take three different values of the power in Eq. Q, 
namely, n = 1, 2 and 3. With the scheme discussed above, we can write the expressions for 
the two-particle twist-3 pion and kaon DAs 



6^ and 
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where P+ = P ■ n. Df denotes the inverse quark propagator, ^ + rrif + Mf{k). The trace 
tr^ denotes the trace over Dirac space. The explicit evaluation of the DAs in Eqs. (fTUI) and 
flll|) is given in Appendix. We set = 140 MeV and mx = 495 MeV for numerical input. 
We assume isospin symmetry with = = 5 MeV. We choose = 150 MeV for the 
strange current quark mass. 



III. NUMERICAL RESULTS 

In the present Section, we provide the numerical results of the twist-3 pion and kaon 
DAs. We fix t he p arameters, using the normalization conditions in Eq. (jH)). The method is 
given m Refs. ^ in detail. 

In Fig. we depict the pseudoscalar-type DAs (0^) for the pion (left) and kaon (right 
For comparison, we also show the asymptotic DA (0Asym.(^) — 1) those of Ref. 



in each panel. Note that the results of Ref. [ll| using the QCD sum rules (QCDSR) were 



ht). 

m 



derived at the renormalization scale /i = 1 GeV, which is rather compatible to our case. As 
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FIG. 1: The results of the pesudoscalar twist-3 pion distribution amphtudes (j)^ (left) and kaon ones 

results of the QCDSR in Ref 

1 is drawn in the long-dashed line 



6^ (right). The dot-dashed curves indicate the results of the QCDSR in Ref. The asymptotic 



one 



Asym. 



expected from isospin symmetry, the pion DAs are all symmetric as shown in the left panel 
of Fig. ^ Examining the dependence of the DAs on the power n in Eq. (j^l), we find a very 
interesting behavior of the DAs with power n. With n = 1 used, the pion DA does not vanish 



at the end points of u and turns out to be very similar to that of Ref. while the pion 
DAs are suppressed in the vicinity of the end points and eventually vanish at the end points 
in the case of n = 2 and 3. In this case, the present results look similar to those of Ref. 0]. 
The reason lies in the fact that for the parameterization of the momentum-dependent quark 
mass with n = 1 there exist nonvanishing terms at the end points after integration over k- 
variables. 

In the right panel, we depict the kaon DAs, (jfj^{u). We see that the kaon DAs are almost 
symmetric in spite of the mass difference between the light and strangeness quarks. Note that 
this behavior was already seen even for the leading twist kaon DA in our previous work . 
It turns out that the present results of the 0^ are very different in shape, compared to 
that of Ref. jl^]. We conclude that, in our model framework, the pion and kaon DAs show 
negligible difference. 

In Fig. El we draw the present results of the pseudotensor twist-3 pion and kaon DAs 
{4>m) whose asymptotic form is ^Xsym. = 6m(1 — u). When n = 1 is used, the 0^ turns 
out to be almost the same as the asymptotic one. However, as n increases, the (j)'^ are 
suppressed rather strongly near the end points, whereas they get narrower and larger in the 
neighborhood of the center {u = 0.5). It can be easily understood from the normalization 



condition for 0^ in Eq. (jS)). Compared to the result of Ref. ll| again, the behavior of the 
present 0^ is of great difference from it. While the present 0^ are humped near the center, 
that of Ref. lUj is rather flat in that region. The pesudotensor twist-3 kaon DAs 0^ look 
similar to the pion ones, as indicated in the right panel of Fig. |2l though they are negligibly 
asymmetric to the right. However, result of Ref. is shifted to the right and slightly 
asymmetric. 

In order to analyze the 0^ and 05vf detail, we compute their Gegenbauer moments. By 
doing that, we can immediately check how much the DAs are deviated from the asymptotic 
ones. In addition, they give information on the strength of flavor SU(3) symmetry breaking 
effects for the kaon DAs. The twist-3 DAs can be expanded in terms of the Gegenbauer 
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FIG. 2: The results of the pesudotensor twist-3 pion distribution amphtudes (p^ (left) and kaon 



ones 0^ (right). The dot-dashed curves indicate the results of the QCDSR in Ref. The 
asymptotic one (pMym = 6m(1 — it) is drawn in the dashed line. 

polynomials: 

oo 

<t>%ii^) = J2<mC'J\0, (12) 

m=0 

oo 

rMiu) = Qu{i -u)Y^ <McTiO, (13) 



m=0 



where ^ = 2m— 1. The orthogonal condition for the Gegenbauer polynomials, the Gegenbauer 
moments for the 0^ and 0^ can be derived as follows: 

<M = (2"^ + 1) duCK'iO (14) 



^0 

Am -I- fi 

where m denotes an order of the Gegenbauer moment. The results of the Gegenbauer 
moments are listed in Table HI First, we consider the pseudoscalar ones, a^. As expected 
from isospin symmetry, it turns out that all odd Gegenbauer moments vanish in the case of 
the pion DAs. Moreover, we see that the Gegenbauer moments are all negative except for 



the case of n = 1. The result of is very similar in general to that of Ref. However, 
the values of 04 and a^.^ are positive in the present work but they are negative in Ref. , 
though in Fig. ^they are seemingly almost the same each other. 

As for the kaon DA, we find that the situation becomes much more interesting. The first 
Gegenbauer moment a^^, which encodes the strength of flavor SU(3) symmetry breaking 
effects, turns out to be all positive in our calculation. Its value is about 0.02. The smallness 
of the Gegenbauer moments explains why the kaon DA becomes almostly symmetric and 
rather flat. Note that the even Gegenbauer moments are much smaller than those of the 
pion DAs. The results of the present work for the kaon DA are much smaller that those of 
the QCDSR 
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We now consider the Gegenbauer moments for the 05^. The overall tendency for the sign 
of the Gegenbauer moments is very similar to that for the 0^. We find that the present 
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results of the are much smaller than those of Ref. p, Lager values of the imply 
that the DA becomes flat, as shown in Fig. |21 Since we obtain all negative 02 if; kaon 
DAs have narrower shapes. Note that all odd kaon Gegenbauer moments are very small 
(~ 10~^). Once again, we see the reduced SU(3) symmetry breaking effects. 
It is also of great interest to study the moments of the DAs defined as follows: 



{rrM= / dur<P%{u), (le) 







(r)^= / ducrMi^). (17) 







The results are hsted in Table m being compared with the QCDSR results from Refs. 0, H, 
HI. We observe that (^^'^'^)^ are comparable with the QCDSR results 0,EI|. We find that 



the values of (C^''^'^)^ are samll in comparison with Refs. (3,lll|, although the even moments 
are still similar. The results of Ref. jsf are slightly larger than the present ones in general. 
As for the moments of the it turns out that the present results are comparable with 



those of Ref. , though the difference between them is getting larger for higher moments. 
The situation is similar when we compare our results with those of Ref. 0. We note that 
the odd moments become all negative. 

Finall y, w e want to discuss the chiral perturbation theory (xPT) Wilson coefficients of 
Refs. p. llC^. based on the present results^. It is of great interest to study them within 
the present model, since the model is chiral symmetric by construction. In Refs. Efll |. 
it was shown that the non-local quark bilinear operators for the twist-3 meson DAs and 
their moments are reconstructed in terms of the chiral field, S = exp (zV2A^"AVFm), to 
the leading order (LO) and next-to-leading order (NLO) as follows: 

(Dm = (OIc^Olo + b^0^i^o\M{P)), (18) 
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1 


0.44296 


0.15081 


0.08421 


0.02033 


0.13715 


0.00170 


0.10748 


-0.00842 


0.05975 


2 


-0.18451 


-0.45274 


-0.26330 


0.02492 


-0.43993 


- 0.02802 


-0.32704 


-0.01548 


-0.16701 


3 


-0.43067 


-0.55590 


-0.17839 


0.02364 


-0.64675 


- 0.03674 


-0.37241 


-0.01999 


-0.09403 


\6] 


0.51578 


0.25454 


0.21624 




0.26310 




- 0.05216 




0.14697 


111] 


0.43726 


-0.07150 


-0.19686 


0.18372 


0.27067 


0.39530 


-0.24693 


0.05497 


-0.24361 


n 




4,7r 
















1 


0.01658 


0.00401 


0.00230 


-0.00359 


-0.00735 


-0.00367 


0.00368 


-0.00280 


0.00188 


2 


-0.09841 


-0.01923 


-0.00368 


-0.00474 


-0.11797 


- 0.00298 


-0.01314 


-0.00068 


-0.00282 


3 


-0.14770 


-0.01596 


-0.00052 


-0.00466 


-0.16072 


- 0.00335 


-0.00949 


-0.00002 


0.00002 


[6] 


0.09400 


-0.00057 


-0.00078 




0.05203 




-0.00048 




-0.00066 


[llj 


0.09787 


-0.00156 


-0.00114 


0.08976 


0.05383 


0.05111 


-0.01501 


0.00734 


-0.00578 



TABLE I: Gegenbauer moments for the twist-3 pion and kaon distribution amplitudes, a' 



The analysis of the xPT Wilson coefficients can be performed by using the derivative expansion within 
the present scheme. The corresponding work will be found elsewhere. 
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where Cm and bm are the xPT Wilson coefficients encoding physics in the scale of A^pt- 
From Eq. p8|) . one can derive the equations concerning the even and odd moments: 





= 0, 








= {ms - 


m)b2m+l,2, 




{enr = 


\? 1 mq = 


,0 + 2rha2m + (2m + m^)/32m, 




{en'/ = 


\? / mq= 


,0 + (m + ms)a2m + (2m + ms)/?2m 


(19) 



where j32m and a2m are the collective parameters of and bm- Note that these equations are 
model- independent, since they are only based on chiral symmetry of the operators [lol |. 
The first equation in Eq. p9|) indicates that the odd moments are all zero for the pion DAs, 
which is true because of isospin symmetry. The second equation represents flavor SU(3) 
symmetry breaking effects in terms of the mass difference between the light and strange 
current quarks, m is defined as (m„ + md)/2 = 5 MeV. As for the even moments, the third 
and fourth equations are given. They can be simplified due to the fact that 0^ and 0^ are 
all zero in the chiral limit as easily verified in Eqs. (llUj) and 



a2n 



{en^ - {en 



K 



'2m 



m + Ills 

(m + m,)(^2m^^ _ 2m(^2m^^ 
mi + ffiiTLs — 2ffi^ 



(20) 
(21) 



It is straightforward to obtain the values of &2m+i by using the results listed in Table ITTl The 
corresponding results are listed in Table IIIIl Since the 62m+i are the expansion coefficients 
for the operators of NLO 0, 3|, we find that the effects from the NLO are getting weaker 
as m increases. In addition to 62™+!, "we also list the calculated values of a2m in Table UTTl 
We verify that the values of /?2m are rather tiny (10~^). Note that /32m are the coefficients 
for the NLO. It turns out that the values of the a2m for the 
they are less than 10~^ for the 



are rather sizeable, while 
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{efn 


{e^K 


{e^K 


{e^K 


{e^K 


{e^K 


1 


0.39240 


0.25445 


0.19072 


0.00678 


0.35162 


0.00416 


0.21840 


0.00292 


2 


0.30873 


0.16742 


0.10820 


0.00831 


0.27468 


0.00338 


0.14142 


0.00160 


3 


0.27591 


0.13666 


0.08164 


0.00788 


0.24710 


0.00263 


0.11663 


0.00081 


[7] 


0.340 ~ 0.359 0.164 ~ 0.211 
















0.52 ±0.03 


0.44 ±0.01 




-0.10 ±0.03 


0.43 ± 0.04 










0.38653 


0.24510 


0.17879 


0.06124 


0.36757 


0.05933 


0.22355 


0.05198 
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{ez 


{en 


{ez 


{Ok 


{e)^, 


{e)^< 


{erK 




1 


0.20568 


0.08992 


-0.00186 


-0.00216 


0.19748 


-0.00162 0.08442 


-0.00131 


2 


0.16626 


0.06122 


-0.00145 


-0.00284 


0.15955 


-0.00179 0.05738 


-0.00123 


3 


0.14936 


0.05030 


-0.00117 


-0.00280 


0.14490 


-0.00184 0.04799 


-0.00125 


[2J 


0.34 ± 0.03 


0.20 ±0.01 




-0.13 ±0.04 0.173 ±0.002 










0.23250 


0.10747 


0.06240 


0.06124 


0.20218 


0.03282 


0.08949 


0.02208 



TABLE II: Moments for the twist-3 pion and kaon distribution amplitudes, (C™)^ • 
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-n nn86n 
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p P p 

04 Ctg 


«2 
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1 


0.26310 0.23258 0.19845 


0.05290 


0.03548 


-0.31503 


2 


0.21968 0.16774 0.12555 


0.04329 


0.02477 


-0.18723 


3 


0.18587 0.12923 0.08890 


0.02877 


0.01490 


-0.14587 



TABLE III: xPT Wilson coefficients 62™+! and a2m defined in Eqs. (UH) and (P7|l . 

IV. SUMMARY AND CONCLUSION 

We investigated the two types of two-particle twist-3 pion and kaon distribution ampli- 
tudes, based on the nonlocal chiral quark model from the instanton vacuum. We considered 
explicitlv flavor SU(3) symmetry breaking effects, employing the modified effective chiral 
action . The current-quark mass dependence was also taken into account for the dy- 

namical quark mass Mf{k) by solving the integral equation for the quark propagator under 
the instanton effects. We employed simple-pole type parameterization for the form factor in 
Mf{k). 

Numerical results of 0^ and (pj^ were given for various cases. We obtained the sym- 
metric pion distribution amplitudes due to isospin symmetry. On the contrary, the kaon 
distribution amplitudes turned out to be asymmetric because of SU(3) symmetry breaking 
effects. However, calculated curves for the kaon DAs were almost symmetric in spite of the 
quark mass differences. 

When power n = 1 for the dynamical quark mass was used, we found that the results 
of the 0^ seem to be very similar to those of Ref. with the same end-point behavior. 
However, as n increased, the results look similar to those of Ref. 0]. As for the the 
results were very similar to that of the asymptotic one, with n = 1 used. As the power 
n increases, the kaon distribution amplitudes were suppressed at the end points and were 
getting narrower. We also investigated the Gegenbauer moments and the moments of the 
distribution amplitudes. 

Finally, we estimated the xPT Wilson coefficients for the twist-3 pseudoscalar meson 
distribution amplitudes, using the present results. These estimations may be useful in ana- 
lyzing the higher-twist meson distribution amplitudes in terms of the chiral field operators. 
A more systematic analysis can be carried out for the xPT Wilson coefficients, the derivative 
expansion being employed. The corresponding analysis will appear soon. 

The investigation on the three-particle twist-3 and twist-4 distribution amplitudes with 
gluon operators is under way. 
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Appendix 



A. Twist-3 pseudoscalar type pseudoscalar meson DA 

The twist-3 pseudoscalar type DA can be written by using Eq. ((T)) as follows: 



U 



V2{P ■ h){mf + 



dr 



-i{2u-l)TP-h 



(0|V'/(rn)^75^9(-™)|A<(P)) 



fng) dr^ -i(2u-i)TP-h 



TT 



{Q\^f{rh)i^,i,g{-rn)\M{P)), 



(22) 



where Fm = /m/V^ stands for the empirical value of the pseudoscalar meson decay con- 
stant. We set rriM to be the mass of the pseudoscalar meson. After integrating over r and 
we arrive at Eq. ()1U|). The trace shown in Eq. (fTUI) can be evaluated in the present 
framework as follows: 



D{k) 
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2 in- 
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2 1 2n' 
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2 1 2" 



nA^ — k'^ 
rUg + Mof{mg 



nA^ 



nA2 -{k- Pf 

2 n 



nA^ 



nA2 - (A; - Pf 



nA^ — k'^ 



2 f2 



X {{k-Pf-ml-2M^f{mg 



nA^ 



nA 



2 H 4n 



nA^ - (A; - Pf 



nA^ — k"^ 

2n 



nA^ 



-\ 4n 



nA2 - (A; - Pf 



X(3k^-kl){ak.~^jf-{(3k.-^g) 



3n 



^VwFai^^k. + ifT{(3k. - 75) 



„ [m/(afc_ - 7/)2" + r]f][mg{l3k_ - 7g)2« + r]g] 



(23) 



where we use the following abbreviations: 

a = P+, (3 = {u-l)P+, 

If = kT + nA^, 'jg = {u- l)m%f + k"^ + nA^, 
6f = k'^ + rnj, 5g = {u - l)m^ + A;^ + m^, 
r^f = Mo/(m;)(nA^)2", r^j = Mof{mg){nA')'\ 

T>f^g in the denominator of Eq. reads: 

'^f,9 = [("^- - ^f,g)i^k^ - 7/,^)^" - '^mj^grif^giak. 
Finally, we have the following expression for the 0^(m): 



(24) 
(25) 



AiNc^/ffJffg{mf + mg)P+ f dk_dk 



77)2 p2 



(27r)3 
10 



(/?fc_-A:j.)(afc_-7/)^"(/3A:_- 7,) 



3n 



-(afc_+7/r(/3A;_-7,r 



(26) 



B. Twist-3 pseudotensor type pseudoscalar meson DA 

One can write the pseudotensor type twist-3 DA, using Eq. Q as follows: 



ruin) 



rrin 



dr 



ml^FM J-oo TT 



X / dve-^'""^^-^ 

10 



{0\^PfiTn)i{pfi - P ■ h)-i^i,g{-Tn)\M{P)). 



(27) 



The trace in the matrix element of Eq. can be derived as in the case of the pseudoscalar- 
type one. Having performed the trace, we finally arrive at 



, ^_ l2iN,{mf + mg)Pl ^ [ dk.d'kr k-jak^ - 7/)'"(/3fc- - 7g)'" 

Jo J "w^. ■ 



(28) 
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